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Decays of the T(5S) resonance to channels with B+ and Z?" mesons are studied using a 23.6 fb ' 
data sample collected on the T(5S) resonance with the Belle detector at the KEKB asymmet- 
ric energy e+e" colHder. The fully reconstructed B+ J/wK^^ ^ J/\j/K*°, B+ D^7t+ 
and D^n^ decays are used to obtain the and Z?" production rates per bb event, 

f{B+) = (67.5 ±3.6 ±4.8)% and f{B^^) = (70.4+^;^ ± 6.2)%. Assuming equal rates to B+ and 
B'^ mesons in all channels produced at the T(5S) energy, we measure the fractions for bb event 
transitions to the two-body and multi-body channels withB+Z" meson pairs, f{BB) = (5.1 ±0.9± 
0.4)%, f{BB*+B*B) = (12.6+}-^ ±1.0)%, f{B*B*) = (34.5; 
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Study ofY(5S) decays to B° and B+ mesons Alexey DRUTSKOY 



New aspects of beauty dynamics can be explored using the large data sample recently col- 
lected by the Belle collaboration at the energy of the T(5S) resonance (also called T(10860)). At 
this energy a bb quark pair can be produced and hadronized in various final states, which can be 
classified as two-body B° channels BJB;, two-body B channels BB, BB*, B*B, 

B*B*, three-body channels BBn, BB*n, B*Bn, B*B*n, and four-body channel BBnn. Here B 
denotes a or a B^ meson and B denotes a or a 5" meson, and the excited states decay to their 
ground states via B* By and B* — > B^y. Fractions and decay parameters for all of these channels 
provide important information about ^-quark dynamics. 

The first study of and B^ mesons at the T(5S) was performed by the CLEO collaboration 
using a 0.42 fb^' data sample, where the fraction of events with meson pairs was found 
to be (58.9 ± 10.0 it 9.2)%. Only two-body B meson channels were observed by CLEO. Several 
theoretical papers are devoted to T(5S) decays to final states with the two-body B*' and B+/° meson 
pairs [0, |3[ ^ IP and with the three-body channels [Q, ^]. The three-body fractions were predicted 
to be about 2-3 orders of magnitude smaller than two-body fractions. Interesting information about 
T(5S) behaviors can be obtained using multi-body decays 

In this analysis we use the data sample of 23.6fb^\ which was taken at the T(5S) center-of- 



mass (CM) energy of ~ 10867 MeV with the Belle detector ||9|] at KEKB [Il0|]. The number of bb 



events in the data sample is Njj^^^'^ = (7.13 ±0.34) x 10^. We fully reconstruct the decay modes 
B+ -^J/\i/K+, B" J/yK*°, B+ D^K+ (D" ^ K+K- and D° K+n+K-K-) and B" ^ D-n+ 
(D K^n^n^). The B decays are reconstructed and identified using two variables: the energy 



difference AE = E^'^ — E^^^ and the beam-energy-constrained mass Mbc = y (£'be^)^ ~ {p'^Y' 
where E™ and are the energy and momentum of the B candidate in the e+e^ CM system, 
and E^jn is the CM beam energy. Two-body and multi-body channels with B+Z" pairs are located 
in kinematically different regions of the Mbc and AE plane, however the central positions for all 
channels are distributed in this plane along a straight line described approximately by the function 
AE = iTiB — Mbc, where ms is the nominal B meson mass. 

The two-dimensional Mbc and AE scatter plots for the B+ J/yK^, J/yK*^, -6+ 
D*'7r+, and B° D modes (with the following //y/ e+e", J /y ^ l^^l^^ ^ K^K^ , 
D" K^n^n^n^, and D K^n^n^ decays) are shown in Fig. 1. Events are clearly concen- 
trated along the B+/" meson production line. 

We used the inclined AE + Mbc — 5.28 projections of the two-dimensional scatter plots for all 
events within the range 5.268 < Mbc < 5.43GeV/c^ to obtain integrated B event yields. To obtain 
the event yields we fit these distributions with a function that includes two terms: a Gaussian to 
describe the signal and a linear function to describe background. Using the fit results, the average 
production rates over studied B modes, /(B+) = (67.5 ± 3.6 ± 4.8)% and /(B°) = (70.4;^^-^ ± 
6.2)%, are obtained. The mean value over B+ and B° modes is /(B+/°) = (68.5+2:9 ± 5.0)%. 
Within errors this rate is in agreement with the CLEO value of /(B+/°) = (58.9 ± 10.0 ± 9.2)%. 

Taking into account the obtained similarity of /(B+) and f{B^) rates, it is natural to assume 
that the B+ and B^ mesons are symmetrically included in all possible channels with B+Z'' meson 
pairs. Therefore the five studied B decays are treated simultaneously everywhere below. First 
the Mbc projections of the two-dimensional Mbc and AE distributions are studied. Only events 
from the signal bands restricted to the ±30 MeV interval in AE (which corresponds to a range of 
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Figure 1: The Mbc and AE scatter plots for the (a) B+ J/vK+, (b) 5° ^ J/\^K*°, (c) B+ 
D^{K+n~)n+, (d) B+ ff' {K+ k+ k-)k+ , and (e) D"7r+ modes. The bands indicate the sig- 

nal regions corresponding to the 5.268 < Mbc < 5.44GeV/c^ and [AS +Mbc - 5.28| < 30MeV intervals. 

(2.5-4.0) c) around the signal linear function (Fig. 1) are used to obtain the signal Mbc distributions. 
To describe combinatorial backgrounds under the signals, we used left and right sideband regions 
corresponding to a shift of 70 MeV in relative to the signal bands. 

Fig. 2a shows the Mbc distributions for the specific two-, three- and four-body channels ob- 
tained from MC simulation for the B'^ D decay. Evidently the two-body channels are 
well separated from each other. The decay matrix elements responsible for the three- and four- 
body decays are not known and the specific three- and four-body channel contributions cannot be 
obtained in a model independent way from the fit of these Mbc distributions. Therefore we re- 
stricted the fit procedure only to the region 5.268 < Mbc <5.348 MeV/c^ to extract the two-body 
channel fractions. To obtain the sum of all three- and four-body channel fractions we used the 
5.348 < Mbc <5.44 MeV/c^ interval for the fit procedure, similar to that used to obtain the full 
/(S+/°) production rate for the entire interval 5.268 < Mbc <5.44 MeV/c^, described above. We 
measured the fractions of bb event transitions to the two-body channels with B+Z*' meson pairs, 
f{BB) = (5.1 ±0.9 ±0.4)%, f{BB*+B*B) = (12.6+};? ± 1.0) %, f{B*B*) = (34.5 +}| ±2.7) %, 
and to the sum of three-body and four-body channels with B+Z" meson pairs, f{B^*''B^*^n{K)) = 
(16.4l}j±1.2)%. The sum of background subtracted Mbc distributions for the five studied B 
decays is shown in Fig. 2b, where the fit results are also superimposed. 

To obtain the three-body channel production rates we additionally used the charged pions 
directly produced in k'^ channels. For all charged pions not contained in the reconstructed 

B candidate, we formed B~/^ 7i:~^ combinations and calculated the values M^^ and AS™^ for the 
missing B meson. To obtain these values we used the energy and momentum of reconstructed Bk 
combination (all values in CM system): £"(5™') = lE^^^^- E{B%) and /^(B™') = p{Bn). 

Figure 2c shows the corrected AE™' ±M^'' - 5.28 projections (m^ = 5.28MeV/c2) for MC 
simulated events where BB 71+ , BB* 7r+ +B*B n^, and B*B* 71+ channels are generated. The recon- 
structed B meson candidates are selected from the signal region within the intervals 5.37 < Mbc < 
5.44GeV/c-^ and lAfi +Mbc — mg| < 30MeV. To improve the resolution we applied additional cor- 
rection requiring the exact equality AE = Mbc — 5.28 for the reconstructed B meson. As we can see 
in Fig. 2c, the BBn^, BB* ;r+ +S*S;r+ and B*B* channel contributions are well separated in 

^mis_^^mis_5 28. 

Finally the AS™** ±M^"^ — 5.28 distribution is obtained in data for the sum of five recon- 
structed B modes (Fig. 2d). We fit the Afi'^' + M^' -5.28 distribution with a function includ- 
ing four terms: three Gaussians with fixed shapes and free normalizations used to describe the 
BBit^, BB* %^ +B*B'ii^, and B*B* 7t^ channel contributions, and a second order polynomial to 
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Figure 2: (a) MC simulated Mbc distributions for the B D decay for (hatched histograms from the 
left to the right) (1) BB, (2) BB*+B*B, (3) B*B* and (4) BBjtTt channels, and also for three-body channels 
BB* n + B*Bn (histogram), BBn (dotted histogram) and B*B* n (dashed histogram). The distributions are 
normalized to unity, (b) Mbc data distribution after background subtraction for the sum of five studied 
B decays (points with error bars) and results of fit (histogram) used to extract two-body channel fractions. 

(c) The A£™'^ + M™** — 5 .28 distribution normalized per reconstructed B meson for the MC simulated — > 
J/\l/K~ decays in the (peaks from the left to the right) BB;r+, BB* ;r+ -|-B*B;r+, and B*B* ;r+ channels. 

(d) The AZi™'** + M™^ — 5.28 data distribution for Z?"/" n'^ combinations for the sum of five studied B modes. 
The curve shows the fit results described in the text. 

describe background. The central positions and widths of the Gaussians are obtained from fits of 
the MC simulated distributions and are fixed in the fit to data. Using the fit results we measured 
the fractions for bb event transitions to the three-body channels with meson pairs, f{BBn) = 
(0.0± 1.1 ±0.2) %, f{BB*n+B*Bn) = (6.8;|^ ±0.7) %, and f{B*B*%) = (LOlj;^ ±0.3) %. The 
three-body rates are calculated assuming the ratio of charged and neutral directly produced pions 
to be 2:1 (Clebsch-Gordan coefficient). It is interesting to note that using the method with the re- 
constructed direct pions, we observe about one half of the rate obtained above for the sum of the 
three-body and four-body channels. This implies a sizable rate for the four-body BBnn channel, 
however no definitive conclusion about this deficit can be reached with the current statistics. 

In conclusion, the production of B+ and mesons is studied at the T(5S) energy. We measured 
the average fractions per bb event for integrated B production, two-body, and multi-body channels. 
The multi-body channels with unexpectedly large fractions are observed for the first time. 
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